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ABSTRACT 

We present simulations of the evolution of a proto-neutron star in which 
kaon-condensed matter might exist, including the effects of finite temperature 
and trapped neutrinos. The phase transition from pure nucleonic matter to the 
kaon condensate phase is described using Gibbs' rules for phase equilibrium, 
which permit the existence of a mixed phase. A general property of neutron 
stars containing kaon condensates, as well as other forms of strangeness, is that 
the maximum mass for cold, neutrino-free matter can be less than the maximum 
mass for matter containing trapped neutrinos or which has a finite entropy. A 
proto-neutron star formed with a baryon mass exceeding that of the maximum 
mass of cold, neutrino-free matter is therefore metastable, that is, it will collapse 
to a black hole at some time during the Kelvin-Helmholtz cooling stage. 

The effects of kaon condensation on metastable stars are dramatic. In these 
cases, the neutrino signal from a hypothetical galactic supernova (distance ~ 8.5 
kpc) will stop suddenly, generally at a level above the background in the SuperK 
and SNO detectors, which have low energy thresholds and backgrounds. This is 
in contrast to the case of a stable star, for which the signal exponentially decays, 
eventually disappearing into the background. We find the lifetimes of kaon- 
condensed metastable stars to be restricted to the range 40-70 s and weakly 
dependent on the proto-neutron star mass, in sharp contrast to the significantly 
larger mass dependence and range (1-100 s) of hyperon-rich metastable stars. 

We find that a unique signature for kaon condensation will be difficult to 
identify. The formation of the kaon condensate is delayed until the final stages of 
the Kelvin-Helmholtz epoch, when the neutrino luminosity is relatively small. In 
stable stars, modulations of the neutrino signal caused by the appearance of the 
condensate will therefore be too small to be clearly distinguished with current 
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detectors, despite the presence of a first order phase transition in the core, fn 
metastable stars, the sudden cessation in the neutrino signal occurs whether it 
is caused by kaon condensation, hyperons or quarks. However, if the lifetime of 
the metastable star is less than about 30 s, we find it is not likely to be due to 
kaon condensation. 

Subject headings: stars: elementary particles - equation of state - stars: interiors 
- stars: neutron 



1. INTRODUCTION 

Proto-neutron stars (PNSs) are formed in the aftermath of gravitational collapse 
supernova, the end state of stars more massive than about 8 M . These objects are 
prodigious emitters of neutrinos of all types which, if detected in terrestrial detectors, could 
reveal details of the supernova mechanism and the properties and composition of dense 
matter. The observation of neutrinos from SN 1987A (Bionta et al. 1987; Hirata et al. 
1987) confirmed the standard scenario for the early evolution of PNSs (Burrows & Lattimer 
1986). At the beginning, this PNS is very hot and lepton-rich and, after a typical time of 
several tens of seconds, the star becomes deleptonized and cold: a neutron star has been 
formed. The detection of neutrinos radiating from the PNS surface is an unequivocal signal 
of the formation of this kind of object, since the direct collapse of the iron core into a black 
hole that might occur, perhaps in relatively extremely massive stars due to accretion, would 
result in a neutrino signal of rather short duration (Burrows 1988). 

One of the chief objectives in modeling PNSs is to determine their internal compositions. 
Many simulations of dense matter predict the appearance of strange matter, in the guise of 
hyperons, a kaon condensate, or quark matter at supernuclear density (Prakash et al. 1997, 
and references therein). An important question is whether or not neutrino observations 
from a supernova could reveal the presence of such matter. One interesting possibility is 
that the existence of strange matter in neutron stars makes a sufficiently massive PNS 
metastable, so that after a delay of 10-100 seconds, the PNS collapses into a black hole 
(Thorsson, Prakash, & Lattimer 1994; Brown 1994; Brown & Bethe 1994; Prakash, Cooke, 
& Lattimer 1995; Glendenning 1995; Keil & Janka 1995; Ellis, Lattimer, & Prakash 1996). 
Such an event might be straightforward to observe as an abrupt cessation of neutrinos 
when the instability is triggered. Once the star becomes unstable, the collapse to a black 
hole proceeds on a time-scale much shorter than the diffusion time scale, and the neutrinos 
still trapped in the inner regions cannot escape (Baumgarte, Shapiro, & Teukolsky 1996; 
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Baumgarte et al. 1996). 

A previous paper (Pons et al. 1999, paper I hereafter) presented calculations of the 
evolution of PNSs, studying the sensitivity of the results to the initial model, the total 
mass, the underlying equation of state (EOS) and the possible presence of hyperons. PNS 
simulations of stars with hyperons, not including hyperon contributions to the opacity, were 
earlier performed by Keil & Janka (1995). Paper I showed that the major effect on the 
neutrino signal before the onset of any possible metastability is the PNS mass: larger masses 
give rise to larger luminosities and generally higher average emitted neutrino energies. In 
addition, it was found that mass windows for metastable models could be as large as 0.3 
M Q , ranging from baryon masses Mb = 1-7 M Q to 2.0 M Q . The lifetimes of hyperonic 
metastable stars decrease with stellar mass and range from a few to longer than 100 s. The 
detection of neutrinos from SN 1987A over a timescale of 10-15 s is thus consistent with 
either the formation of a stable PNS or a metastable PNS containing hyperons, as long as 
its mass was less than about 0.1 M below the maximum mass for cold, catalyzed hyperonic 
matter. Larger PNS masses would lead to a collapse to a black hole on a timescale shorter 
than that observed. 

A similar situation could be encountered if the EOS allowed the presence of other 
forms of "exotic" matter, manifested in the form of a Bose condensate (of pions or kaons) 
or quarks (Prakash, Cooke, & Lattimer 1995; Ellis, Lattimer, & Prakash 1996; Steiner, 
Prakash, & Lattimer 2000). As suggested some years ago (Kaplan & Nelson 1986), kaons 
obtain an effective mass which decreases with density in dense matter. As a consequence, 
the ground state of hadronic matter at high density might contain a kaon condensate (Brown 
et al. 1992; Muto & Tatsumi 1992; Thorsson, Prakash, & Lattimer 1994; Maruyama et al. 
1994), which would result in a much softer EOS (Tatsumi & Yasuhira 1999; Glendenning 
& Schaffner-Bielich 1998; Pons et al. 2000, hereafter paper II). However, a high lepton 
fraction, which exists in the early PNS evolution, suppresses the formation of kaons just as 
hyperon formation or the appearance of quarks is impeded (Prakash, Cooke, & Lattimer 
1995; Prakash et al. 1997). The appearance of kaons after deleptonization produces a 
softening of the EOS, which can destabilize a sufficiently massive star, causing a collapse 
into a black hole. 

In this paper we concentrate on models obtained in the context of kaon condensation 
and explore their minimal and maximal effects. One important difference between models 
containing hyperons and those containing a kaon condensate stems from the first order 
phase transition that often exists in the latter case. For some models of kaon-nucleon 
interactions, mixed phase regions can be present. Previous works (Glendenning 1992; 
Prakash, Cooke, & Lattimer 1995; Glendenning & Schaffner-Bielich 1999, paper II) have 
emphasized that satisfying the Gibbs' phase rules for equilibrium in situations when more 
than one conserved charge (in our case, baryon number and electric charge, lepton number 
being constrained by beta equilibrium) is present can radically change the pressure-density 
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relation compared to the case in which a Maxwell construction is employed. Knorren, 
Prakash, & Ellis (1995); Schaffher & Mishustin (1995) have emphasized that the threshold 
density for the appearance of strange particles, hyperons and kaons, are sensitive to the 
rather poorly known interactions in dense matter. In models in which hyperons appear at 
lower densities than kaons, the role of kaons can be relatively small. On the other hand, 
should it turn out that the interactions of the charged hyperons are relatively repulsive, 
only neutral hyperons may be present and kaons can play a dominant role. 

Although the possibility of PNS metastability was first discovered in the context of 
kaon condensation in neutron star matter (Thorsson, Prakash, & Lattimer 1994), a full 
simulation of a PNS evolution with a kaon condensate has not been performed so far. 
Baumgarte, Shapiro, & Teukolsky (1996) have previously studied the hydrodynamical 
collapse of metastable PNSs containing kaon condensates. They were, however, mostly 
interested in the dynamical collapse to a black hole, not in the PNS evolution. They 
employed the kaon EOS from Thorsson, Prakash, & Lattimer (1994) in which the 
transition from normal matter to the phase containing kaons was described by a Maxwell 
construction. They also rescaled the neutrino mean free paths by a factor of a thousand 
to accelerate the quasi-static part (i.e., the Kelvin-Helmholtz phase) of the evolution, in 
order to save computation time. This produces little effect on the collapse to a black hole, 
since neutrinos are essentially trapped in the short time of collapse. However, their results 
for the Kelvin-Helmholtz stage can only be considered as rough estimates, because of this 
over-simplification. 

In this paper, we perform the first evolutionary simulations of PNSs with kaon 
condensates during the Kelvin-Helmholtz epoch by using Henyey-like evolution and 
flux-limited diffusion modules as described in paper I. As in paper I, we ignore the effects of 
accretion, which are generally believed to be small after a successful supernova explosion, 
expected to occur within 1 s of core bounce. We use a field-theoretical EOS that includes 
kaon condensation (paper II) using boson exchange interactions and employ Gibbs' phase 
rules for multicomponent matter in which baryon number, electric charge, and lepton 
number are conserved during the evolution. We use the opacities developed by Reddy, 
Prakash, & Lattimer (1998); Reddy et al. (1999) and Burrows & Sawyer (1998, 1999) 
in the nucleon sector and account for the modifications that arise from the presence of 
kaon-condensed matter (Thorsson et al. 1995; Ji & Min 1998; Reddy 1998). 

Our presentation is organized as follows. §2 provides an outline of the basic features of 
the EOS of matter containing a kaon condensate. In §3 we describe the thermodynamical 
properties of kaon-condensed matter at finite temperature and with varying amounts of 
leptons likely to be encountered in the evolution of a PNS. §4 contains a description of the 
neutrino opacities in kaon condensed matter. Results from numerical simulations of the 
Kelvin-Helmholtz phase of the evolution of a PNS with a kaon condensate are presented in 
§5. Here we also explore the sensitivity of the results to the initial mass and consequences 
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associated with metastable PNSs which collapse into a black hole. The detectability of 
black hole formation in the Super-Kamiokande (SuperK) and SNO detectors from possible 
future galactic supernovae are investigated in §6. Our conclusions are given in §7. 



2. EQUATION OF STATE OF KAON CONDENSED MATTER 

The EOS of kaon-condensed matter including the effects of trapped neutrinos and finite 
temperature for several different field-theoretical models for both the nucleon-nucleon and 
kaon-nucleon interactions was studied in paper II. Here we summarize the basic relations 
required for the computation of the EOS. Since our purpose in this paper is to explore 
whether a discriminating neutrino signal would be observed from kaon condensed matter, we 
present results for cases in which a first order phase transition occurs at zero temperature. 

To model the baryonic phase, we use a field-theoretical description in which baryons 
interact via the exchange of a, u, and p mesons. Including leptons, the Lagrangian density 
is given by (Serot & Walecka 1986), 

L = Lh + Li 

= Y,B(-i^d, - gvB-fwp - gpBl^^ ■ t - M*)B 

B 

i 

Here, B are the Dirac spinors for baryons, t is the isospin operator, and the gs and ms are 
the couplings and masses of the mesons. The sums include baryons B = n,p, and leptons, 
I — e~, fx~ and v { . The field strength tensors for the u> and p mesons are W^ u = d fM u u — d v ui^ 
and B M j, = d^h^ — d u h p , respectively. The potential U(a) represents the self-interactions of 
the scalar field and is taken to be of the form (Boguta & Bodmer 1977) 

1 i 

U (a) = -bm B (g aN a) + -c(g aN a) , (2) 

where m# is the bare nucleon mass. 

In the mean field approximation the thermodynamic potential per unit volume is 
(Serot & Walecka 1986) 

^ = \mW + U(a) - \mlul - \m% - 2T £ J || In (l + e^*—)) . (3) 
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Here the inverse temperature is denoted by (3 = l/T and E* = \/k 2 + M* 2 with 

M* = m B — g a cr denoting the nucleon effective mass. The chemical potentials are given by 

Hp = v v + gjjj Q + \9pbo ; fl n = v n + 9^0 ~ \g P b • (4) 
Using Qjj-, the thermodynamic quantities can be obtained in the standard way. 

The leptons are included in the model as noninteracting particles, since their 
interactions give negligible contributions compared to those of their free Fermi gas parts. 

For the kaon sector, we use a Lagrangian that contains the usual kinetic energy and 
mass terms along with the meson interactions (Glendenning & Schaffner-Bielich 1999, paper 
II). Kaons are coupled to the meson fields through minimal coupling; specifically, 

L K = V^K+WK- -m%K + K- , (5) 
where the vector fields are coupled via the standard form 

£> M = d M + ig U KUn + WpKl^p ■ t (6) 
and m* K = tuk — \g a K^ is the effective kaon mass. 

In the mean field approach, the thermodynamic potential per unit volume in the kaon 
sector is (paper II) 



oo 



^ = Um 2 ^K - {}* + *o) 2 ) + Tj -*1L [l„(l _ e -««-i0) + i n( i _ e -^))] , (7) 

o 

where X = g^K^o + gpxbo, the Bose occupation probability = (e^ x — l) -1 , 

w ± = + m *z ±x , f = 93 MeV is the pion decay constant and the condensate 
amplitude, 9, can be found by extremization of the partition function. This yields the 
solution 9 = (no condensate) or, if a condensate exists, the equation 

m* K = fi K + X (8) 

where \ik is the kaon chemical potential. In beta-stable stellar matter the conditions of 
charge neutrality 

^2 ^ bHb - n e - n K = (9) 



and chemical equilibrium 



to = biVn ~ Qiito ~ Hv e ) (10) 

= Hn-Hp (11) 
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are also fulfilled. 

The kaon condensate is assumed to appear by forming a mixed phase with the baryons 
satisfying Gibbs' rules for phase equilibrium (Gibbs 1876). Matter in this mixed phase is in 
mechanical, thermal and chemical equilibrium, so that 

V l = V U , T 1 — T u , f4 = l4 J , (12) 

where the superscripts I and II denote the nucleon and kaon condensate phases, respectively. 
The conditions of global charge neutrality and baryon number conservation are imposed 
through the relations 

xq* + {i-x)q" = o 

+ (1 - xWb = n B , (13) 

where \ denotes the volume fraction of nucleonic phase, q the charge density, and % 
the baryon density. We ignore the fact that the phase with the smallest volume fraction 
forms finite-size droplets; in general, this would tend to decrease the extent of the mixed 
phase region. Further general consequences of imposing Gibbs' rules in a multicomponent 
system are that the pressure varies continuously with density in the mixed phase and that 
the charge densities must have opposite signs in the two phases to satisfy global charge 
neutrality. We note, however, that not all choices of nucleon-nucleon and kaon-nucleon 
interactions permit the Gibbs' rules to be satisfied (for an example of such an exception, 
see paper II). The models chosen in this work do allow the Gibbs' rules to be fulfilled at 
zero and finite temperatures and in the presence of trapped neutrinos. 

The nucleon-meson couplings are determined by adjusting them to reproduce the 
properties of equilibrium nucleonic matter at T = 0. We use the numerical values used by 
Glendenning & Moszkowski (1991), i.e., equilibrium density n = 0.153 fm -3 , equilibrium 
energy per particle of symmetric nuclear E/A = —16.3 MeV, effective mass M* = 0.78M, 
compression modulus Kq = 240 MeV, and symmetry energy a sym = 32.5 MeV. These values 
yield the coupling constants ga/m a = 3.1507 fm, g w jm w = 2.1954 fm, g p /m p = 2.1888, b = 
0.008659, and c= -0.002421. 

The kaon-meson couplings g a x and g^K are related to the magnitude of the kaon 
optical potential Uk at the saturation density n of isospin symmetric nuclear matter: 

U K (n ) = -g aK o-(n ) - sw^oK)- (14) 

Fits to kaonic atom data have yielded values in the range —(50 — 200) MeV (Friedman, 
Gal, & Batty 1994; Friedman et al. 1999; Waas & Weise 1997; Ramos & Oset 2000; Baca, 
Garcia- Recio, & Nieves 2000). We use g^ic = fl'wAf/3 and g p x = g P N '/2 on the basis of simple 
quark and isospin counting. Given the uncertainty in the magnitude of \Uk\, consequences 
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for several values of \Uk\ were explored in paper II. Moderate values of \Uk\ generally 
produce a second order phase transition and, therefore, lead to moderate effects on the 
gross properties of stellar structure. Values in excess of 100 MeV were found necessary for 
a first order phase transition to occur; in this case kaon condensation occurs at a relatively 
low density with an extended mixed phase region, which leads to more pronounced effects 
on the structure due to a significant softening of the EOS. 

3. COMPOSITION AND STRUCTURE OF PNS WITH KAON 

CONDENSATES 

In Figure 1 we show the variation of pressure with baryon density for three different 
choices of entropy per baryon, s, and lepton fraction Y L = (n e + n M + n u )/nb. The 
optical potential Uk is set to —120 MeV. In an evolving PNS, s and Yl are not constant, 
but these three choices are reasonable approximations to the ambient conditions before 
(s = 1, Yl = 0.35) and after (s = 2, Y u = 0) the deleptonization stage, as well as 
for a cold, neutrino-free (s = 0, Y v = 0) neutron star. The dotted line corresponds to 
lepton-rich matter, the dashed line corresponds to hot neutrino-free matter and the solid 
line corresponds to cold neutrino-free matter. Results obtained from Gibbs' construction 
are shown by thick lines while those for the pure phases are shown by thin lines. 

Applying Gibbs' rules eliminates the region of negative compressibility (this is more 
evident in the case of cold matter), and also decreases the critical density for the onset of 
kaon condensation. Above the threshold density for condensation, the kaon concentration 
builds up rapidly, which results in a significant softening of the EOS. The effects of trapped 
neutrinos are similar to those obtained for hyperons in Prakash et al. (1997) and paper 

I. The most notable feature due to the presence of trapped neutrinos is that the critical 
density for kaon condensation is much higher in the case neutrinos are trapped than in the 
case of neutrinos having left the star. Notice that the extent of a Maxwell construction 
would also be largely reduced in neutrino-trapped matter, and the differences between both 
kinds of mixed phase, Maxwell and Gibbs, are very small. 

The phase boundaries of the different phases are displayed in Figure 2 in a Yl-ub 
plane for an optical potential Uk of -100 MeV (left) and -120 MeV (right), respectively. 
The nucleonic phase, the pure kaon matter phase, and the mixed phase are labelled I, 

II, and III, respectively. Solid lines mark the phase transition at zero temperature and 
dashed lines mark the phase transition at an entropy per baryon of s = 1. Note that finite 
entropy effects are small and do not affect significantly the phase transition density. The 
dash-dotted line shows the electron fraction Y e as a function of density in cold, catalyzed 
matter (for which Yl = Y e ), which is the final evolutionary state. The region to the left 
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of this line corresponds to negative neutrino chemical potentials and cannot be reached 
during normal evolutions. The solid and dashed lines, which separate the pure phases from 
the mixed phase, vary roughly linearly with the lepton fraction. Also notice the large, 
and nearly constant, densities of the boundary between the mixed phase III and the pure 
kaon phase II. These densities, for the cases shown, lie above the central densities of the 
maximum mass stars, so that region II does not generally exist in proto-neutron stars (see 
Figure 5 below). The effect of increasing the lepton number is to reduce the size of the 
mixed phase (which in fact shrinks to become a second order phase transition for Yl > 0.4 
and Uk = —100 MeV) and to shift the critical density to higher densities. A similar effect 
is produced by decreasing the magnitude of the optical potential. 

Figure 3 shows the baryon (upper panel) and lepton (lower panel) chemical potentials 
for the three ambient conditions employed in Figure 1. In neutrino-free matter, we observe 
a rapid decrease of the electron chemical potential produced by the gradual substitution 
of electrons by kaons. When neutrinos are trapped, /3-equilibrium prevents the electron 
chemical potential from decreasing, because the neutrino chemical potential increases with 
density. In Figure 4 the compositions are for the same ambient conditions of Yl and s as 
in Figure 1. A remarkable effect is the very high proton fraction reached at high densities, 
higher than the neutron fraction at densities above 0.8 fm -3 . The figure also displays how 
both thermal and lepton-trapping effects displace the critical density to higher values, 
the effects of having neutrinos in the system being more important than those of finite 
temperature. 

In Figure 5 we show the baryon mass (M B , upper panels) and gravitational mass (M G , 
lower panels) versus central baryon density, for Uk = —120 MeV (left panels) and -100 
MeV (right panels), respectively. The onset of kaon condensation is shown by the diamond 
on each curve. In each case, the maximum mass of the cold neutrino-free star is below 
the maximum mass of the hot star, irrespective of whether it has neutrino trapping. This 
figure indicates the existence of a mass window for metastability (1.7 < Mb < 2.0 M 
for U K = -120 MeV and 2.0 < M B < 2.2M Q for U K = -100 MeV) in which a PNS will 
become unstable at the end of its Kelvin-Helmholtz stage. A similar effect was found in 
paper I due to the presence of hyperons in matter. It should be emphasized that this figure 
uses constant s and Yl profiles within the star. In fact, these quantities are not constant in 
the interior, so the precise mass window depends upon the initial profile and the evolution. 
We have found, however, that the mass window exhibited in this figure is a reasonable 
approximation. 

The main effect of lowering \Uk\ is a reduction in the size of the metastability window 
and its displacement to higher values, and this is shown in Figure 6. The value of Uk must 
be limited to values greater than about -126 MeV in order that the binary pulsar mass 
constraint, 1.44 M , not be violated. For values of Uk greater than about -80 MeV, it 
is also clear that metastable models cannot exist. Another interesting feature is that the 
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maximum masses of hot, neutrino-free models are nearly the same (for XJ K = —120 MeV) 
or even larger (for Uk = —100 MeV) than those of lepton-rich configurations, preventing 
metastable stars from collapsing upon deleptonization and delaying the onset of instability 
until the end of the cooling stage. This means that kaon-induced instabilities cannot occur 
very quickly, as can be the case with hyperon induced collapses. 



4. NEUTRINO OPACITIES 

For neutral and charged current processes involving nucleons and leptons, we utilize 
the recent advances made in the works of Reddy, Prakash, & Lattimer (1998); Reddy et 
al. (1999); Burrows & Sawyer (1998, 1999). In these works, neutrino cross sections for 
arbitrary degeneracies were calculated, including the effects of interactions and collective 
excitations in beta stable isospin asymmetric matter. These opacities form the basis of our 
baseline calculations against which modifications to the opacities due to the presence of a 
kaon condensate are assessed. 

In bulk equilibrium and in the degenerate limit, the factors by which the normal 
reaction rates must be modified in the presence of a charged kaon condensate are (Thorsson 
et al. 1995): 

v e + n(K) ->• n(K) + e~ : ^ sin 2 9 tan 2 9 C 
v e + p(K) -* p(K) + e" : sin 2 9 tan 2 9 C 

u e + n(K) -> p(K) + e~ : cos 2 (#/2) (15) 

Here, n(K) denotes an excitation which is a superposition of a neutron and a proton, and 
reduces to a free neutron or proton in the absence of a charged kaon (or pion) condensate, 
9 denotes the amplitude of the condensate and 9c is the Cabibbo angle (sin# c = 0.23). 
The factors in this equation are all less than unity. The opacity for the neutral current 
scattering processes involving nucleons in the presence of a kaon condensate has also been 
investigated recently by Ji & Min (1998) and by Reddy (1998), who included the SU(3) 
singlet contribution. The general result is that a significant reduction in the cross section is 
caused by the presence of a kaon condensate formed in a bulk medium. 

However, the situation is radically different for the opacity due to a kaon condensation 
in a mixed phase, in which finite-size effects may be important (Glendenning & 
Schaffner-Bielich 1998; Reddy, Bertsch, & Prakash 2000; Christiansen, Glendenning, & 
Schaffner-Bielich 2000). Reddy, Bertsch, & Prakash (2000) found that coherent scattering 
of neutrinos off kaon-condensed droplets in the mixed phase increases the neutrino cross 
section, relative to bulk nucleonic matter, by factors of 10-20. The overall enhancement, 
however, is sensitive to the droplet size, which depends on the surface tension. 
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We consider two extreme models: one in which the cross section or opacity is 
significantly smaller (20 times) than that of the neutrino-nucleon processes to account for 
a condensate in bulk matter, and one in which the cross section is significantly larger (20 
times), to account for possible effects of a droplet phase. We apply these factors to the 
overall opacity, while in fact the reduction factor for the bulk condensate should apply only 
to the volume occupied by the kaon-condensed phase. The mixed phase occupies only a 
small fraction (~ 0.2) of the total volume even in the most massive PNSs, which mitigates 
considerably the large reductions in the opacity in this case. We wanted to ensure that our 
models will unmistakably bracket the true situation. Our results conclusively demonstrate 
that even large modifications to the opacity caused by kaon condensation do not have a 
pronounced effect on neutrino signals. 



5. EVOLUTION 

We now proceed to the evolution of a PNS in which a kaon condensate forms. We will 
focus on the case in which Uk = —120 MeV, for which the effects of kaon condensation 
are relatively pronounced. Calculations (Prakash et al. 1996; Muto, Tatsumi, & Iwamoto 
2000) of the growth rate for a kaon condensate indicate a timescale of less or about 10~ 4 
s at the temperatures of interest (above a few MeV). Since this is much shorter than 
the Kelvin- Helmholtz time-scale (several s), we will assume matter to reach chemical 
equilibrium instantaneously. However, for hydrodynamical studies, such as those involving 
collapse to a black hole or pulsations due to core quakes, the kaon condensation formation 
timescale might have to be considered. 

Initial entropy per baryon s and lepton fraction Yl profiles are taken to be the same 
as in paper I, where we studied the evolution of PNSs containing either pure nucleons or 
nucleon/hyperon mixtures. The numerical code used to perform the simulations was also 
described in paper I; the code is based on a Henyey-like scheme in which the structure and 
transport equations are solved in sequential steps and then corrector steps are taken until 
convergence is reached. 

For M B < 1-6 M , the central density does not exceed the critical value for kaon 
condensation (for the set of parameters employed), and the evolution is identical to that of 
pure nucleon (np) models described in paper I. For stars with 1.6 < M#/M < 2.05, a kaon 
condensate will form during the evolution. We can classify the PNSs with a kaon condensate 
core in two subclasses: i) stable, those with Mb < 1.75 M Q , which is the maximum mass 
of cold, neutrino-free neutron stars (see Figure 5); ii) metastable, those configurations with 
1.75 < M B /M & < 2.05. 

The evolution of stars containing kaon condensates [npK) deviates from those of np 
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stars when the mixed phase with a condensate forms at the core. The size of this mixed 
phase grows on a diffusion time scale, i.e., the average time it takes for the neutrinos 
produced during kaonization to escape the inner core. In Figure 7, we show profiles of 
particle fractions as a function of radius for a stable star with Mb = 1.7M Q towards the 
end of the Kelvin- Helmholtz stage (t = 60 s). At this time, the star is at the beginning of 
its long-term cooling phase during which its final catalyzed state is achieved. Notice that 
the size of the core is about 3 km, but it contains only a small fraction (?» O.O5M ) of the 
star's total mass. 

In Figure 8, the evolution of stable (1.7 M ) and metastable stars are compared. Both 
the central baryon number densities and kaon fractions Yk are displayed. In each case, 
the time at which kaon condensation occurs is indicated by a diamond. Asterisks mark 
the times at which the evolution of metastable stars could not be further followed in our 
simulations, i.e., when a configuration in hydrostatic equilibrium could not be found. At 
this time, the PNS is unstable to gravitational collapse into a black hole. For the stable 
star, kaons appear after about 40 s. Thereafter, the star's central density increases in a 
short interval, about 5 s, until a new stationary state with a mixed phase is reached. The 
evolution of the metastable stars is qualitatively different, inasmuch as the central density 
increases monotonically from the time the condensate appears to the time of gravitational 
collapse. 

It is interesting that the lifetimes in all cases shown lie in the narrow range 40-70 s 
(see Fig. 8). They decrease mildly with increasing M B . Further insight is obtained by 
examining Figures 9 and 10. In Figure 9 the evolution of the central baryon density and 
kaon fraction are shown as functions of the central value of Yj,. The initial contraction 
of the star occurs at nearly constant lepton fraction, as can be seen from the vertical 
trajectories of the central density. This is followed by a stage in which the lepton fraction 
decreases while little variation in the central density occurs, until kaon condensation takes 
place (marked by diamonds in the figures). Thereafter, both the central density and the 
kaon fraction increase until the PNS reaches its final stage, either the cold-deleptonized 
configuration or the unstable configuration that collapses to a black hole. 

Although the densities and lepton fractions for which condensation occurs are different 
for stars with different masses, the typical time scales are similar, as is evident from Figure 
10 in which the evolution of the electron and neutrino concentrations, Y e (top panel), Y u 
(middle panel), respectively, and the temperature (bottom panel) at the star's center are 
shown. For all cases, the lepton concentrations exhibit a plateau stage until about 30-40 
s, and remain high enough to prevent the formation of a kaon condensate. Once the 
condensate is formed, Y v may increase as kaons gradually replace electrons in the inner core. 

Notice that the plateau stage is also present in the evolution of stars containing only 
nucleons, as was shown in paper I. Until a kaon condensate appears, the evolution of stars 
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considered in this work proceeds similarly to that of a nucleons-only star because the 
number of thermal kaons present is rather small. The plateau is caused by an increase in 
neutrino opacity due to high temperatures, which reduces the leakage of neutrinos from 
the core, in turn maintaining a high lepton fraction. Consequently, kaonization of matter 
is delayed until the end of this epoch, due to the dependence of the threshold density 
on Y L . This explains why kaon-condensed stars of different masses become unstable at 
approximately the same time. This is qualitatively different from metastable hyperon-rich 
stars studied in paper I, for which the lifetime is a much more sensitive function of mass. 
This is displayed in Figure 11, in which the lifetimes as a function of Mb for stars containing 
hyperons (npH) and npK stars are compared. In both cases, the larger the mass, the 
shorter the lifetime. For kaon-rich PNSs, however, the collapse is delayed until the final 
stage of the Kelvin-Helmholtz epoch, while this is not necessarily the case for hyperon-rich 
stars. 

Strangeness appearing in the form of a mixed phase of strange quark matter also leads 
to metastability. Although quark matter is also suppressed by trapped neutrinos (Prakash, 
Cooke, & Lattimer 1995; Steiner, Prakash, & Lattimer 2000), the transition to quark 
matter can occur at lower densities than the most optimistic kaon case, and the dependence 
of the threshold density on Yl is less steep than that for kaons. Thus, it is an expectation 
that metastability due to the appearance of quarks, as for the case of hyperons, might be 
able to occur relatively quickly. Calculations of PNS evolution with a mixed phase of quark 
matter, including the possible effects of quark matter superfluidity (Carter & Reddy 2000) 
are currently in progress and will be reported separately. 

An interesting question concerns how short the lifetime of a kaon-condensed metastable 
star can be. For the case Uk = —120 MeV, this time is constrained by the fact that 
metastability disappears for Mb ~ 2.05M Q . For smaller magnitudes of the optical potential, 
the lifetimes are larger since the critical densities are larger; these densities are reached only 
after longer times. Conversely, increasing the magnitude of Uk decreases the lifetime. The 
magnitude of the optical potential is limited, however, by the binary pulsar mass constraint 
that the maximum mass must exceed 1.44 M , or \U K \ < 126 MeV for our interaction. In 
this case, the lifetime is found to be about 40 s. This leads to the interesting conclusion 
that should metastability from a PNS with a timescale of, say, 30 s or less be observed, it 
would be evidence for the existence of hyperonic or quark matter rather than for a kaon 
condensate. 
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6. LUMINOSITIES AND SIGNALS IN TERRESTRIAL DETECTORS 

The time dependence of the neutrino signal from a supernova is the only observable 
means through which the physical processes occurring in the high density inner core of a 
PNS may be inferred. Our study is more applicable to times longer than approximately 
1 s after core bounce, after which effects of dynamics and accretion become unimportant. 
Studies of the neutrino signal during the first second, during which approximately 1/3 of 
the neutrino energy is emitted, require more accurate techniques for neutrino transport 
(Messer et al. 1998; Yamada, Janka & Suzuki 1999; Burrows et al. 2000) coupled with 
recent developments in supernova simulations (Mezzacappa et al. 2000; Rampp & Janka 
2000). 

For simplicity, we limit our study to the signal produced by electron antineutrinos 
absorbed onto protons in a pure water detector, for which the cross section is o(E u ) = a El 
where E v is the v e energy in MeV and a = 9.3 x 10~ 44 cm 2 . This was the major 
reaction observed during the SN 1987A neutrino burst and it is by far the process with the 
highest rate in the newer detectors. Assuming that the neutrinos leaving the PNS have a 
Fermi-Dirac spectrum with zero chemical potential, knowledge of the neutrino luminosities 
and average temperatures permit the signal from a supernova to be calculated according to 
(paper I) 

~dt ~ A^ M F 3 (0) T " Lu ■ (16) 

In this equation, L D is the antineutrino luminosity (assumed to be 1/6 of the total neutrino 
luminosity), n p = 6.7 x 10 28 free protons per kiloton of water, D is the distance to the 
supernova, M. is the fiducial mass of the detector, and is the redshifted neutrino 
temperature. Also, -F 3 (0) = 77r 4 /120 is an ordinary Fermi integral and Gi(E t h,T) denotes a 
modified, truncated Fermi integral: 

POO 

G t (E th ,T)= / dzz l W(zT)(l + e z )- 1 (17) 

with E th being the detector threshold and W(E) the detector efficiency We studied the 
signal observed in KII and 1MB, whose masses, thresholds and efficiencies are taken to be 
the same as in Lattimer & Yahil (1989), except that 6 ktons is used for the 1MB mass, as 
discussed by (Schramm & Brown 1990). We have assumed thresholds and fiducial masses 
of 5.6 (2) MeV and 22.4 (1.7) ktons for the SuperK (SNO) detectors, respectively, and also 
assumed a detector efficiency of unity above the threshold for both. Although we assumed 
that PNS neutrinos have a zero chemical potential spectrum, the formula (16) does not 
depend very sensitively on this assumption Lattimer & Yahil (1989). 



The luminosity Lbg corresponding to a given background rate below which the signal 
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is lost, (dN/dt)BG, can be obtained from inverting equation (16): 



fdN\ D/lOkpc 2 
Ibg = 105.0 — \' n P ( F T P °° , 18 
V / BG (M/kton) v " ' ' v ; 

where F(T?°) also depends on the threshold and efficiency of each detector. We have 
obtained the following analytical fits 

1 + exp [2.0(1.2 — T)] 



F(T) = < 



25.5(T-0.75) 

(-5.27T + 2.57T 2 )- 1 

1 + exp [2.25(1.2 - T)} 
23.3T 

1 + exp [6.0(0.4 — T)] 



KII, T > 0.75 
1MB, T > 2.05 
SuperK 
SNO 
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23.3T 

where T is in MeV. In general, we will assume that the limiting background rate is about 
0.2 Hz (one count every 5 s), since the corresponding time between counts is then about 
10% of the entire signal time. A better estimate would require an in-depth statistical 
analysis that is beyond the scope of this paper. 

In Figure 12 the evolution of the total neutrino energy luminosity is shown for different 
models. Notice that the drop in the luminosity for the stable star (solid line), associated 
with the end of the Kelvin-Helmholtz epoch, occurs at approximately the same time as for 
the metastable stars with somewhat higher masses. In all cases, the total luminosity at the 
end of the simulations is below 10 51 erg/s. The two upper shaded bands correspond to SN 
1987A detection limits with KII and 1MB, and the lower bands correspond to detection 
limits in SNO and SuperK for a future galactic supernova at a distance of 8.5 kpc. The 
width of the bands represents the uncertainty in the v e average energy because of the 
limitations of the flux-limited diffusion approximation, as discussed in paper I. The times 
when these limits intersect the model luminosities indicate the approximate times at which 
the count rate drops below the background rate (dN/dt)BG = 0.2 Hz. 

The poor statistics in the case of SN 1987A precluded a precise estimate of the PNS 
mass. Nevertheless, had a collapse to a black hole occurred in this case, it must have 
happened after the detection of neutrinos ended. Thus the SN 1987A signal is compatible 
with a late kaonization-induced collapse, as well as a collapse due to hyperonization or to 
the formation of a quark core. More information would be extracted from the detection of 
a galactic SN with the new generation of neutrino detectors. 
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In SNO, about 400 counts are expected for electron antineutrinos from a supernova 
located at 8.5 kpc. The statistics would therefore be improved significantly compared to 
the observations of SN 1987A. A sufficiently massive PNS with a kaon condensate becomes 
metastable, and the neutrino signal terminates, before the signal decreases below the 
assumed background. In SuperK, however, up to 6000 events are expected for the same 
conditions (because of the larger fiducial mass) and the effects of metastability due to 
condensate formation in lower mass stars would be observable. 

It is interesting to assess the role of neutrino opacities in the presence of kaon 
condensates. For this purpose, we calculated the luminosities with two extreme assumptions 
concerning the opacities, namely, the opacity is (i) 20 times smaller than the baseline 
opacity, corresponding to the bulk case (Thorsson et al. 1995; Ji & Min 1998), and (ii) 
20 times larger than the baseline opacity, corresponding to finite-size effects in the mixed 
phase (Reddy, Bertsch, & Prakash 2000). 

In Fig 13, we show the temporal evolution of the central baryon density (top panel) 
and the total luminosity (bottom panel), for models with and without the correction to 
the opacity. In addition, we also show results corresponding to the case of pure np matter, 
in which the kaon condensate is not allowed. The magnitude of the opacity in the mixed 
phase chiefly controls the growth rate of the central density (top panel), and hence of the 
condensate, with larger opacities causing longer delays. The effect of reduced opacities on 
the luminosities (bottom panel) is barely distinguishable from the baseline case. In the case 
of enhanced opacities, for the stable star, the luminosity is first reduced by about 10-20 % 
for 10 s following the appearance of kaons at t ~ 45 s, but eventually exceeds that of the 
baseline case. Nevertheless there is only a small change in the lifetime produced by even 
large opacity changes. The metastable star in Figure 13 with 2M Q collapses before much 
information about the physics in the core reaches the neutrinosphere. It is also remarkable 
that the largest differences in the luminosities between the models with and without a 
condensate, but for the same opacity, are smaller than about 2%. 

It is important to note, however, that by the time significant effects due to a kaon 
condensate become visible (t > 50 s), our treatment of neutrino transport becomes suspect 
since neutrinos in the interior are entering the semi-transparent regime. This is reflected in 
the location of the neutrinosphere, which begins to fall well below the surface of the star 
for these times, and in the total optical depth of the star, which is no longer large. These 
affect our estimates of both the neutrino signal and the effective detector background. 

The expected neutrino signals in the SuperK detector from a supernova at 8.5 kpc are 
compared in Fig. 14 for a sample of stable PNSs with different compositions and masses. To 
show the effects of compositional changes, stars with Mb = 1.7 M for np, npH and npK 
cases are shown, and to show the effects of small changes in the stellar mass, an Mb = 1.8 
M for the case of npK is also displayed. The upper panel shows the count rate and the 
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lower panel shows the integrated number of counts. It is clear that, for the same stellar 
mass and for stable models, the effect of having a kaon condensate core is smaller than the 
effect of hyperons: the np case is indistinguishable from the npK case. A change in the 
assumed mass of the PNSs of about 0.1 M has about twice the effect on the number of 
counts as the change produced by introducing hyperons. Therefore, unless an independent 
method of identifying the PNS mass is available, small differences in the assumed PNS mass 
can mask compositional effects. 

Furthermore, although a clear signature of the presence of strangeness in PNSs occurs 
due to the collapse of metastable models, there is no unique feature that establishes the 
presence of a kaon condensate. If the lifetime of the metastable star is less than about 30 s, 
however, it appears unlikely that the instability is due to kaon condensation. 

The main conclusions arising from these comparisons is that the presence of a kaon 
condensate would be difficult to establish in an unambiguous fashion. For stars that are 
stable, the small modulations in the time dependence of the neutrino luminosity or count 
rate due to kaon condensation can be easily masked by differences in composition (i.e., 
the presence of hyperons) or a small difference in the assumed PNS mass compared to a 
nucleons-only case. For stars that are metastable, the sudden cessation in the neutrino 
signal is similar for instabilities whether caused by kaons, hyperons or quarks. 



7. SUMMARY AND CONCLUSIONS 

We have studied the effect of kaon condensation on PNS evolution, by employing an 
EOS which includes the effects of finite temperature and neutrino trapping. The phase 
transition from pure nucleonic matter to the kaon condensate phase is described by means 
of Gibbs' rules for phase equilibrium which permit a mixed phase. In the models explored 
here, the central densities are not large enough to allow a pure condensed phase to exist. 

We can classify stars of different masses in three main groups: i) stars in which the 
central density does not exceed the critical value for kaon condensation, ii) stars that can 
form a mixed phase core at the end of the Kelvin-Helmholtz epoch but remain stable, and, 
iii) kaon condensed metastable stars that become unstable and collapse to a black hole at 
the end of the Kelvin-Helmholtz epoch. For stars in which the effects of kaon condensation 
are small (this could be either because the star's mass is low enough to permit only a 
very small region of the mixed phase or because condensation occurs as a weak second 
order phase transition), the differences of the predicted neutrino signal compared to PNSs 
composed of pure nucleonic matter are very small. 

Our calculations show that the variations in the neutrino light curves caused by the 



-19- 



appearance of a kaon condensate in a stable star are small, and are apparently insensitive 
to large variations in the opacities assumed for them. Relative to a star containing only 
nucleons, the expected signal differs by an amount that is easily masked by an assumed 
PNS mass difference of 0.01 — 0.02 M . This is in spite of the fact that, in some cases, a 
first order phase transition appears at the star's center. The manifestations of this phase 
transition are minimized because of the long neutrino diffusion times in the star's core and 
the Gibbs' character of the transition. Both act in tandem to prevent either a "core-quake" 
or a secondary neutrino burst from occurring during the Kelvin-Helmholtz epoch. 

Observable signals of kaon condensation occur only in the case of metastable stars 
that collapse to a black hole. In this case, the neutrino signal for a star closer than 
about 10 kpc is expected to suddenly stop at a level well above that of the background 
in a sufficiently massive detector with a low energy threshold such as SuperK. This is in 
contrast to the signal for a normal star of similar mass for which the signal continues to 
fall until it is obscured by the background. The lifetime of kaon-condensed metastable 
stars has a relatively small range, of order 50-70 s for the models studied here, which is in 
sharp contrast to the case of hyperon-rich metastable stars for which a significantly larger 
variation in the lifetime (a few to over 100 s) was found. This feature of kaon condensation 
suggests that stars that destabilize rapidly cannot do so because of kaons. 

We determined the minimum lifetime for metastable stars with kaons to be about 
40 s, by examining the most favorable case for kaon condensation, which is obtained by 
maximizing the magnitude of the optical potential. The maximum optical potential is 
limited by the binary pulsar mass constraint, which limits the star's maximum gravitational 
mass to a minimum value of 1.44 M . Therefore, should the neutrino signal from a future 
supernova abruptly terminate sooner than 40 s after the birth of the PNS, it would be more 
consistent with a hyperon- or quark-induced instability than one due to kaon condensation. 

It is important to note that the collapse to a black hole in the case of kaon condensation 
is delayed until the final stages of the Kelvin-Helmholtz epoch, due to the large neutrino 
diffusion time in the inner core. Consequently, to distinguish between stable and metastable 
kaon-rich stars through observations of a cessation of a neutrino signal from a galactic 
supernovae is only possible using sufficiently massive neutrino detectors with low energy 
thresholds and low backgrounds, such as the current SNO and SuperK, and future planned 
detectors. 

This work was supported in part by the Spanish DGCYT grant PB97-1432 and the 
U.S. Department of Energy under contract numbers DOE/DE-FG02-87ER-40317 (JAP 
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Fig. I.- 

Total pressure as a function of baryon density for different lepton fractions and 
entropy per baryon in matter with a kaon condensate. Thin lines show results from the 
pure phases while thick lines are results obtained employing Gibbs' phase rules. 

Fig. 2.- 

Phase boundaries between pure nucleonic matter (I), pure kaon condensed matter 
(II) and a mixed phase (III) in the Yl-tib plane for Uk = —100 MeV (left panel) and 
Uk = —120 MeV (right panel). The solid line corresponds to s = and the dashed line to 
s — 1. The dashed-dotted line shows the baryon density as a function of the lepton fraction 
for s = 0, neutrino-free (Yj, = Y e ) matter. 

Fig. 3.- 

Baryon (top panel) and lepton (bottom panel) chemical potentials for the physical 
conditions shown in Figure 1. 

Fig. 4.- 

Particle concentrations as a function of baryon density % for the physical conditions 
shown in Figure 1. 

Fig. 5.- 

Baryon (top) and gravitational (bottom) mass versus central baryon density for 
Uk = —120 MeV (left) and -100 MeV (right), respectively, for the physical conditions 
shown in Figure 1. Diamonds indicate the density above which kaon condensation occurs. 

Fig. 6.- 

The maximum baryon mass for the physical conditions shown in Figure 1 as a function 
of the optical potential Uk- Metastable stars cannot exist when the cold, catalyzed 
maximum mass exceeds that of neutrino trapped matter. The triangle indicates the binary 
pulsar gravitational mass constraint Mq°- x = 1.44 M . 

Fig. 7.- 

Particle concentrations as a function of the radius of a PNS containing a kaon 
condensate core. The figure is a snapshot at t = 60 s, near the end of the Kelvin-Helmholtz 
stage. Notice the increase of proton concentration in the presence of kaons. 
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Fig. 8.- 

Evolution of the central baryon density % and the kaon/baryon fraction Y K for stars 
with different baryonic masses. Solid lines correspond to a stable star with M B = 1.7M Q ; 
stars with larger masses are metastable. Diamonds indicate when kaon condensation occurs 
at the stellar center, and asterisks denote when metastable stars become gravitationally 
unstable. 

Fig. 9- 

Evolutionary tracks of the central densities and kaon fractions as functions of Yl for 
stars of various baryon masses. Diamonds and asterisks have the same meaning as in Figure 
8. 

Fig. 10.- 

Evolution of the central temperature (bottom) and neutrino (middle) and electron 
(top) fractions for stars of various baryon masses. Diamonds and asterisks have the same 
meaning as in Figure 8. 

Fig. 11.- 

Lifetimes of metastable stars as a function of the stellar baryon mass. Solid lines show 
results for PNSs containing kaon-condensates and dashed lines show the results of paper I 
for PNSs containing hyperons. 

Fig. 12.— 

The evolution of the total neutrino luminosity for stars of various baryon masses. Line 
and symbol designations are the same as in Figure 8. Shaded bands illustrate the limiting 
luminosities corresponding to a count rate of 0.2 Hz in all detectors, assuming a supernova 
distance of 50 kpc for 1MB and Kamioka, and 8.5 kpc for SNO and SuperK. The width of 
the shaded regions represents uncertainties in the average neutrino energy from the use of a 
diffusion scheme for neutrino transport. 

Fig. 13.- 

The evolution of the central baryon density (top) and total neutrino luminosity 
(bottom) for two stars with M B = 1.7 and 2.0 M , for different assumptions concerning 
composition and opacity. Solid lines refer to the baseline model of a kaon-condensed star 
which assumes opacities corresponding to pure nucleonic matter. Dotted lines refer to stars 
containing only nucleons. Dashed (dot-dashed) lines refer to models of kaon-condensed 
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stars in which the opacity is increased (decreased) by a factor 20. Asterisks indicate when 
metastable models become unstable to gravitational collapse. The shaded bands are as in 
Figure 12. 

Fig. 14.— 

The expected neutrino signal in the SuperK detector from a supernova at 8.5 kpc for 
PNSs with different compositions and masses. For the case Mb = 1.7 M Q , np, npK and 
npH compositions are compared. For the npK composition, stars with two masses of (1.7 
and 1.8 M Q ) are compared. The upper panel shows the count rate and the lower panel 
shows the integrated number of counts. 
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